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Abstract 

We describe a technique for estimating the floating body po- 
tentials of partially-depleted silicon-on-insulator (PD-SOI) 
circuits under steady switching activity and under initial 
activity after a long period of quiescence. The approach is 
based on a unique state diagram abstraction of the PD-SOI 
FET that captures all of the essential device physics. This 
picture yields a simple analytic model of the body voltage 
which is used within the context of a prototype transistor- 
level static timing analysis engine. Results are presented 
that demonstrate the accuracy of the analytic body-voltage 
model and the reduction in delay uncertainty possible with 
this technique. 

1 Introduction 

Silicon-on-insulator (SOI) technology has long found niche 
applications for radiation-hardened or high-voltage integrated 
circuits. Recently, SOI has emerged as a technology for high- 
performance, low- power deep-sub micron digital integrated 
circuits [1, 2], For digital applications, fully-depleted devices 
have been largely abandoned in favor of partially-depleted 
technology because of the difficulty in controlling the thresh- 
old voltage of fully-depleted thin-film transistors. Partially- 
depleted SOI (PD-SOI) has two main advantages for digital 
applications: the reduction of parasitic source-drain deple- 
tion capacitances and the reduction of the reverse- body ef- 
fect in stack structures and pass- transistor logic. 

At the device and circuit level, however, the floating 
body effect in partially-depleted SOI (PD-SOI) poses "un- 
certainties" in the body potential, and hence the threshold 
voltage of a FET, while the parasitic bipolar effect can result 
in noise failures[5]. Design margining required to offset un- 
certainties in body voltage can erode potential performance 
advantages under nominal operation. For circuit styles in 
which noise margin is strongly determined by threshold volt- 
age (e.g. dynamic circuits), considerable overdesign for noise 
can also result from conservative body- voltage margining. 
Previous circuit-level modelling work on PD-SOI has fo- 
cussed on device issues [3, 4] or delay and noise effects due 
to the floating-body effect evident for particular circuits 
under periodic stimulus [5, 6] (pulse stretching, frequency- 
dependent delay time). In this paper, we present the first 
techniques to quantify floating- body effects over tens of mil- 
lions of transistors through static analysis. 

We consider characterizing the body- voltage uncertainty 
of PD-SOI devices using knowledge of switching activity. 
Each FET of a circuit is analyzed to determine the minimum 
and maximum possible body voltage (the body- voltage un- 
certainty) under both "initial-condition" and "steady-state" 



operation. Initial-condition operation holds for circuits that 
are quiescent for a long time before undergoing switching ac- 
tivity. Steady-state body- voltage estimation takes into ac- 
count average switching behavior of the circuit that has been 
consistently present for a long time. Two techniques are con- 
sidered for steady-state body- voltage estimation. The first 
assumes only that the circuit is under reasonably steady 
switching activity while the second technique, which has 
less uncertainty, requires knowledge of signal probabilities 
and arrival times. This body- volt age estimation technique 
is applied to transistor-level static timing analysis [7], but 
can be extended to static noise analysis [8] It can also be 
used in circuit simulation to avoid the need for long simula- 
tion runs to establish steady-state body- voltage values. 

In this paper, we work with BSIM3SOI [9] models for an 
IBM partially-depleted SOI technology described elsewhere 
[10]. Devices have a 0.25fira effective channel length, 5 — 
nm gate oxide, 350 — nm back oxide, and 140 - nm thin 
silicon film. Two supply voltage are considered - 1.0V and 
2,5V. While the detailed results presented here apply to this 
technology, the techniques are generally applicable to any 
PD-SOI technology. 

In Section 2, a state-diagram abstraction of the PD- 
SOI FET which simplifies the device physics determining 
the body voltage is presented. In Section 3, a simple an- 
alytic model is derived which can be used to accurately 
predict the body voltages in PD-SOI circuits under both 
initial-condition and steady-state operation. Section 4 de- 
scribes a prototype transistor-level static timing analysis 
engine which incorporates these body voltage characteriza- 
tions. Some timing analysis results are presented in Section 
5. Section 6 concludes and offers direction for future work. 

2 State diagram view of body interactions 

The body potential of a PD-SOI FET is determined by ca- 
pacitive coupling of the body to the gate, source, and drain, 
by diode currents at the source- body and drain- body junc- 
tions (including gate-induced drain leakage (GIDL) [11]), 
and by impact ionization currents produced by current flow 
through the device (sometimes referred to as the on-state 
impact ionization current). To model the switching his- 
tory determining the body voltage of a particular device, 
we use the state diagram abstraction shown in Figure 1 
(for an nFET). The states denoted with solid circles rep- 
resent "static" states, valid logic conditions for a FET at 
the end of a cycle. For example, state 1 corresponds to 
the case in which the gate is high and both the source and 
drain axe low. States 6a and 6b, denoted with dashed cir- 
cles, are "dynamic" states since for CMOS designs that do 
not draw dc current, these states will be present only tran- 
siently during switching events. Arrows indicate possible 
state transitions produced by switching events in the cir- 
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cuits containing these FETs. These switching events can 
represent transitions from the logic state at the end of the 
previous cycle to the logic state at the end of the current 
cycle or can represent hazards that occur transiently within 
a cycle. The states 1, 2, 5a, and 5b are distinguished in Fig- 
ure 1 with dashed boxes because transitions between these 
states propagate timing delays. These are also the state in 
which FETs would be susceptible to noise (glitches). States 
5a and 5b can usually be treated equivalently as state 5; sim- 
ilarly states 6a and 6b can be treated equivalently as state 
6. The state diagram of the pFET is the "dual" of Figure 
1, in which the gate is high rather than low in states 3, 4, 
and 5; and low rather than high in states 1, 2, and 6. 



state diagram abstraction of Figure 1. We are implicitly as- 
suming that the digital circuits under consideration are part 
of a discrete-time system characterized by a cycle (usually 
defined by the action of a clock) with a cycle time t cycle - 
This enables us to consider the behavior of the body over a 
long period of time to be determined by an "average" cycle 
repeated over and over. Such an average cycle is shown in 
Figure 2, divided into a series of time slices which char- 
acterize the amount of time per cycle on the average that 
the FET spends in state i. Of course, the time slices sum to 
the cycle time: 



(1) 




Figure 1: State diagram for a PD-SOI nFET. 

In analyzing the device physics determining the body 
potential, it is convenient to distinguish "slow" and "fast" 
processes. Fast processes can change the body potential on 
time scales on the order of or less than the cycle time, while 
slow processes require time scales much longer than the cycle 
time (up to milliseconds) to affect the body voltage. There 
are two fast mechanisms at work: switching transitions on 
the gate, source, or drain which are capacitively coupled 
to the body, and forward- bias diode currents across source- 
body and drain-body junctions with voltages exceeding the 
diode turn-on voltage. The coupling "kicks" associated with 
the first mechanism occur for each transition in Figure 1 and 
are completely reversible on "fast" time scales; that is, if one 
begins in state 1 and traverses the state diagram, returning 
to state 1 on a time scale comparable or less than the cycle 
time, the body voltage on return will be the same as the ini- 
tial body voltage (a simple statement of charge neutrality). 
The reversibility of the kicks is, of course, dependent on not 
triggering the fast irreversible discharge of the body through 
a forward-biased junction. In practice, this discharge mech- 
anism is only triggered on activation of a circuit after a long 
period of quiescence. Once the discharge happens (usually 
within the first few cycles of operation), it is not triggered 
again as long as the circuit is under steady switching ac- 
tivity. This discharge, when associated with a transition to 
state 5, can produce parasitic bipolar leakage. 

The slow processes involve charging or discharging the 
body through reverse-biased or very weakly forward-biased 
diode junctions and through impact ionization. These leak- 
age currents give each state a comparatively slow (much 
longer than the cycle time) relaxation to a target dc body 
voltage for each state in Figure 1. 



3.1 Body-voltage model 

In characterizing the body voltage, we choose a reference 
state of state 2 for the nFET and state 1 for the pFET. 
The discontinuities between the slices occur as a result of 
the capacitive coupling kicks and are characterized by the 
voltage differences, or displacements, di, (see Figure 2) be- 
tween the body voltage in state i and the body voltage in the 
reference state (the reference body voltage, which we denote 
as Vg**). The particular choice of reference states we have 
made ensures that all of the di are negative for the nFET 
and positive for the pFET. In addition, each state is relax- 
ing to the target dc values, The time constants for this 
relaxation are denoted as n. Except for the fast discharge 
associated with a source-body or drain- body junction that 
becomes strongly forward biased, these time constants are 
much larger than f cycle, and any voltage change during a 
single time slice would be imperceptible in Figure 2. From 
this simple picture, one can relate the body voltage at the 
end of the cycle v n+1 to the body voltage at the beginning 
of the cycle v n by: 

w n+1 = » n e'=> +2 -( (s.— *)(!-« ' /ri )e ,=,+1 

1 = 1 

(2) 

The steady-state solution of this difference equation (in the 
approximation that the n are much greater than f cycle) is 
given by: 
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3 Determining the body potential 

In this section, we introduce a simple analytic model for 
the floating body potential of a PD-SOI FET based on the 



Figure 2: An average cycle for an nFET which when re- 
peated over and over models the behavior of the body over 
a long period of time. 



532 



It is important to note that both the displacements d, 
and the relaxation times r; are body* bias dependent; that 
is, both di and n are functions of Vj e/ so that Equation 3 
must be solved self-consistently. The body- volt age depen- 
dencies of di and r, are modelled as simple polynomial curve 
fits to circuit simulation results, since the relationships are 
too complex to motivate physically-based formulae. (We 
actually fit to the logarithm of t\, since r; varies over sev- 
eral orders of magnitude.) We find that both the n and d, 
are independent of device width, W. di is independent of 
W because the dominant components of capacitance scale 
proportionately with W. Moreover, n is independent of W 
because both the body currents and body capacitance scale 
proportionately with W, 

We let Vj' tro denote the steady-state reference voltage 
in state j, V/ tro = sj - dj. Also, above (below) a certain 
value of Vg ! for the nFET (pFET), one or both of the 
source-body or drain-body junctions become strong forward 
biased (which we consider to be at biases greater than 0.6 
V). We refer to the reference voltage at which this occurs as 

yfvmard for j 

We note that, at a supply voltage of 2.5V in our example 
technology, V r / oru,ard = V 6 ' cro because of the dominating 
effect of the on-state impact ionization current at this supply 
voltage. This current is so large that a strongly forward- 
biased junction is required to balance it in steady-state. This 
means that whenever the device spends any time in state 6 
at this supply, the steady-state reference body voltage is 
immediately (within a microsecond) pinned at yJ orward - 

K 6 

The model of Figure 1 and Equation 3 provides for both 
initial-condition and steady-state body-voltage analysis. We 
consider four types of estimation. As is traditionally done 
in static analysis, each has associated with it minimum and 
maximum possible values that bound all of the remaining 
unknowns, the body voltage uncertainty. 

Full-uncertainty body-voltage estimation. This ap- 
plies when there is no knowledge of the switching activity 
of the circuit and one must be assured that the uncertainty 
covers all possible stimulus and history. One considers a 
given state ; to be accessible if i*** > 0 for that state. 
The dynamic state 6 is considered accessible if more than 
one of the states 1, 2, and 5 are accessible, since switch- 
ing events between these states almost always involves tran- 
siently passing through state 6. Let A denote the set of all 
accessible states, and let A at atic denote the set of all acces- 
sible static states. In the full-uncertainty case, one assumes 
that the FET can spend an indeterminant amount of time 
in any accessible state. In this case, the maximum (min-r 
imum) body voltage is the maximum (minimum) value of 
Vf* ro across all accessible states. For the nFET and pFET, 
the minimum and maximum body voltages are then given by 
WUn = min,^ Vr° and (V^) max = max j€ , V/"°, 
respectively. 

Initial-condition body- voltage estimation. This ap- 
plies for times in an initial period of switching activity af- 
ter a long period of quiescence. In this case, one assumes 
that the FET has reached a dc steady state in any ac- 
cessible static state (i. e., states 1 through 5). For the 
nFET and pFET, the minimum and maximum body volt- 
ages are then given by (V£ e/ ) min = min i6 ^, (a| . c Vf cro and 
(V£ c/ )m«* = maXj €>i , (a(io Vf m , respectively. We do not 
consider state 6 here because it can only be reached as part 
of a switching event and does not contribute to the body 
voltage "near* dc quiescence. The uncertainty comes about 



because we have assumed that we do not know the specific 
quiescent state of the circuit. 

Accessibility steady-state body-voltage estimation. 
This applies to the case that the circuit is under unknown 
but steady switching activity. In this case, we take ad- 
vantage of the fact that for the nFET (pFET), KI e/ can 
not be greater (smaller) than the largest (smallest) value 
of W oru,or<l in any accessible state j\ otherwise, the body 
would rapidly discharge. This, of course, assumes that all 
accessible states are visited with reasonable frequency. Fur- 
thermore, for the nFET (pFET), the smallest (largest) V£ e/ 
can be is the smallest (largest) value of Vf ero in any acces- 
sible state j. Therefore, in this analysis, if we know the 
FET is under reasonably steady switching activity in which 
all accessible states are visited, the minimum and maxi- 
mum body voltages for the nFET are given by (V^ e, ) m in = 
mm jeA Vf< r ° and (^) WM = min ;€M v/ oru ' ord , respec- 
tively. For the pFET, the minimum and maximum body 
voltage are given by {V?') min = max i6 ^ V/°™ ard and 
(V£ e/ )mo* = max> € >i Vf* ro t respectively. The uncertainty 
of accessibility estimation can be reduced if one is further 
willing to restrict the allowable waveforms to those meeting 
known timing requirements and known signal probabilities. 

Detailed steady-state body- voltage estimation. This 
applies to the case in which the circuit is under switching 
activity that can be quantified logically with known signal 
probabilities and temporally with known arrival times. We 
tighten up the steady-state body- voltage uncertainty pro- 
vided by accessibility analysis by determining the maximum 
and minimum possible V£ ef values from Equation 3 with 
consideration of all of the i** f . To maximize V£ e/ , one 
must maximize the time in those states with higher Vf* ro ; 
similarly, to minimize V£ e/ , one must maximize the time in 
those states with lower Vf cro . To formalize this, we define 
the rank of the state i as an integer rf indicating the priority 
of states to maximize the body voltage. We can also define 
a complementary rank r» for minimizing V£*'; these prior- 
ities are simply the reverse of the maximum case. One's 
ability to favor the of given states in detailed body- 
voltage estimation is limited by additional constraints de- 
rived from static timing analysis and stochastic techniques, 
as described later in this section. 

3.2 Determining the accessible states and allowable val- 
ues of t i ss 

To determine the accessible states A for a given FET as 
well as to determine allowable values of t \ u for detailed 
steady-state analysis, we need information about the circuit 
environment of each transistor, both logical and temporal. 

We characterize the logical environment of each FET by 
a set of signal probabilities which determine the possible 
states of the source, gate, and drain of the transistor at the 
end of a cycle. For the nFET, we define the probabilities (or 
conditional probabilities) P{G), P{D\Z!) t P(D|G), P(S\G) } 
P(S|G), P(S\G) = P{D\G), and P(S\G) = P(D\G). P(S\V) t 
for example, is the probability that at the end of a cycle 
the source is driven high given that the gate isjow. There 
are analogous probabilities for the pFET with G taking the 
place of G. Accessibility of a state can be immediately de- 
termined from these probabilities. For example, nFET state 
1 is accessible if {P(G) > 0) /\(P(S\G) < 1). 

To calculate the t \ n necessary for detailed steady-state 
analysis, we need to use these FET signal probabilities to 
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calculate the transition probabilities, Pi—>, the probability 
that a given FET in a given cycle is transitioning from state 
* to state j. To find these transition probabilities, we be- 
gin by noting that the cycle-to-cycle transitions through the 
state diagram in Figure 1 represents a Markov process with 
six-by-six transition matrix A. Assigning different transi- 
tion matrices to the minimum and maximum cases yields 
the difference equations (P t mas )k+i = A max (/ > / nttX ) fc and 
(Pr in )k+i = A m <„(/r <n )*, where the matrix A ma ^ min 
has each column given by 

/max /min max f min ^max / min _max / min max/min max/min\T 
(Pi P2 Pz Pa P5a Phb ) 

(Pi)k is the probability of being in state t at the end of 
cycle h. pi is the probability of making a transition to 
state t and follow directly from the source, gate, and drain 
signal probabilities. Forjexample, for the nFET, p mo * = 
(1 -P(G))P{S\G)P(D\G) while p? in = (1 - P(G))(l - 
P(S\G))(l - P(D\G)). The maximum (minimum) case as- 
sumes that the floating node condition on the source or drain 
takes a high (low) voltage value. Diagonalizing A (triv- 
ially) and finding the eigenvector associated with eigenvalue 
1 (normalized so that the sum of the elements of the vec- 
tor is 1) gives the steady-state values of the Pi. From these 
probabilities, one can calculate a set of thirty-six transi- 
tion probabilities for both the minimum or maximum cases: 

pmaxfmin pmax fm\n max/min 

^> ~~ i P} 

Calculating the t\ n for detailed steady-state analysis 
involves determining not only the probability of making a 
transition in a given cycle but the fraction of the cycle time 
(*i~* fc ) tnat be spent in each state t as part of a given 
transition j -* k. In particular, we seek the values of tj~"* 
to maximize and minimize the body voltage among the set 
of possible waveforms. To do this, we require temporal in- 
formation, in particular early and late arrival times (rising 
and falling) for the source, gate, and drain of the FET un- 
der consideration (the target FET). We denote these arrival 
times for the early case as: 

• S;t;l V >D' r ul V ,G) a a n y - EaxUest times the source, drain, 
gate of the FET can be driven high. 

• S/ol/ V ' ^/ait"' ^rili y * Earliest times the source, drain, 
gate of the FET can be driven low. 

There are comparable arrival times associated with the late 
state: S l r VA, sfiK. D W> *nd <%&l For the 

source and dram arrival time, we are assuming that the tar- 
get FET is off. The details of how these are obtained in the 
context of static timing analysis is described in Section 4. 
Each transition has associated with it a set of arrival times 
necessary to make that transition, a transition set, denoted 
as Ti^j. For example, for the 1 -* 2 transition for the nFET, 
the associated arrival time set is T\~i = {S r i 8e , D r i$c}* 
For 1 4 for the nFET, the associated arrival time set 
is = {GfaUiSrim Drise}. We can define max and min 
operators which act on the transition sets. max(7i— j) re- 
turns the largest of the early arrival times in the transition 
set, while min(7i— j) returns the smallest of the late arrival 
times in the transition set. 

To indicate the states of Figure 1 involved in a cycle 
and to handle the possibility of hazards, we can denote the 
waveform in a cycle (in this case involving a transition from 

i to j) using the transition notation t ™* it ^ n j. In this 
cycle, a hazard to state k occurs as part of the transition. 
The transition notation must involve only static states and 



indicates the amount of time spent in each of these static 
states as part of the transition. Specifically for this exam- 
ple, <p = max(7U fc ). tp J = min(7Uj) - max(7^ fc )> 
and ty = t cycle - min(7i-*j). Hazards are introduced 
when they act to increase (in the case that we are seeking 
the maximum body voltage) or decrease (in the case that 
we are seeking the minimum body voltage) the steady-state 
body voltage that would result from the particular wave- 
form being repeated indefinitely. Specifically, a state k can 
be inserted between states t and j according to one of the 
following cases. 

Case I. If (r k > r;) /\{r k > r>) /\{mzx{Ti-. k ) < 
min(7i— j)), the state k can be inserted between i and j as a 
hazard to increase the steady-state body voltage. The cycle 

after this insertion is i ^5-* k ^ n j. In this case, increas- 
ing the amount of time in state k at the expense of states 
i and j increases the body voltage. If (f k > fi)A(^ fc ^ 
fj) A(max(7i— *) < min (7*— j)), the state k can be inserted 
between i and j as a hazard to decrease the body voltage. 
The cycle after this insertion is the same as in the maximum 
case. However, in this case, increasing the amount of time 
in state k at the expense of states i and j decreases the body 
voltage. 

Case II. If (r fc > rj)/\{r k < n) f\{mm{Ti~ k ) > 
min(75— >)) A(min(7i— >) > min(7i— *)), the state k can be 
inserted between t and j as a hazard to increase the body. 

voltage. In this case, the initial cycle is t j. After 

insertion, it is t ^V* k j. State k is inserted only if 
adding it does not decrease the time in state i. If (r k > 
rj)/\fjk < ri) A(min(7^ fc ) > min^j)) /\{mm{T k ^) > 
min(7i— *)), the state k can be inserted between t and j as a 
hazard to decrease the body voltage. The cycle before and 
after insertion is the same as in the maximum case. 

Case III. If {r k > n) f\(r k < r > )/\(max(T fc _ J ) < 
max(7i^>)) y\(max(7it-j) > max(Ti^ k )) } state k can be 
inserted between i and j as a hazard to increase the body 
voltage. Initially the cycle is t ^i* j. After insertion, it 
is i k ™ x j. State k is inserted only if adding it 
does not decrease the time in state /. If (r> > r,) /\(r k < 

r>) A( mu (7*-v') < max(7;-.j)) A(max(^fc->) > max(Ti~ k 
state k can be inserted between i and j as a hazard to de- 
crease the body voltage. The cycle before and after insertion 
is the same as in the maximum case. 

These preliminaries lead to a straightforward algorithm 
for determining the t k ~* J to maximize or minimize the body 
voltage for a given transition. For the maximum case (the 
minimum case is the same except the complementary ranks 
fi are used): 

1. If rj > r k} then the starting cycle is j f ^i n k else the 
8 tar ting cycle is j k, 

2. Find the accessible state k' with the largest rank (dif- 
ferent from j and k and not previously inserted or at- 
tempted) that can be inserted between states j and k, 
If no such state exists, then exit with the current cycle. 
However, if such a state exists then the new cycle is 

. max . . min , . max . ; mo* . . min . » min , 

either j — + k —+ k t j -*■ k — ► ft, or j Jb, 
depending on which case led to the insertion. 

3. Repeat step 2 for each transition in the current aver- 
age cycle, Repeat this until no further refinement is 
possible. 
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For the purposes of detailed steady-state body- voltage es- 
timation, we assume that if state 6 is accessible, exactly 
one switching event through 6 occurs per cycle and assign a 
small switching time to each . We find that, in practice, 
the detailed results are not very sensitive to the exact value 
of t l <T J * J U6t tnat there is some nonzero time in state 6. 

From the determined above, we can calculate an 
effective amount of time (<*") on the average per cycle that 
the FET is in the state t: 

^effynaxfmin ^ ^ pmoi/m<n^j->fcjinaj/min 

>.* 

Obviously, there are two sets of values, one to maximize 
and one to minimize the body voltage. These determine 
the uncertainty of the detailed steady-state body voltage 
estimate. 

4 Static timing analysis 

We have incorporated this body voltage characterization 
into a prototype static transistor-level timing analysis en- 
gine, SOISTA [for SOI Static Timing Analyzer]. The design 
is partitioned into channel-connected components (CCCs) 
for analysis, as is traditionally done in static transistor-level 
tools [7]. SOISTA utilizes a breadth-first search (BFS) of the 
resulting timing graph and propagates signal probabilities 
using assumptions of spatial and temporal independence, 
borrowing from similar techniques in static power analysis 
[12]. Signal probabilities can be easily propagated from CCC 
inputs to outputs using basic probability theory [13] or BDD 
analysis [12, 14]. Uncertainty is built into the timing analy- 
sis to account for variations in the effective gate input capac- 
itance depending on the switched state of the gate (loading 
uncertainty). 

Once the signal probabilities and arrival times are known 
at the inputs of a CCC, these probabilities are translated 
into the FET signal probabilities and arrival time values 
for detailed body- volt age estimation. We first consider the 
calculation of the signal probabilities. If we let i and j de- 
note two channel nodes of the CCC, then similar to [15], 
we can define the kth path P*j as one connection of FETs 
between i and j. We can also define a path function f P * 

as a Boolean function indicating whether the kth path is 
conducting. Let n, denote a controlling nFET gate input 
function in the path, and let pi denote a controlling pFET 
gate input function in the path. Then, the path function is 
given by fp* = A^ei* ni A Ap j€ i* - ^ If there are N 
paths between t and >, then the total path function fp itj is 
given by fp.j = \Jy P h Jp*.* The path probability P(Pi,j), 

the probability that at the end of the cycle > the path from i 
to j is conducting, follows from elementary probability the- 
ory or BDD analysis [12]. The source and drain conditional 
probabilities required for the detailed body voltage estima- 
tion are given by specific path probabilities. For example, 
P(D\G) is the path probability between the drain node and 
Vdd with the gate of the target transistor low. 

We next consider calculating the FET arrival times which 
determine the temporal circuit environment of each FET. 
G /"i <y » G r"i V i G /°an> G YaU are determined directly from 
the CCC arrival times. To determine S'?* v , D%£ y , D^l 
$rVAt we trace all paths from the channel node (node t) to 



Vdd (node 1) with the target transistor ofL 
VA.i : S*°;l v = min (Vn ilW € : (GJSX) 
VP..1 : Sj?i: = max (Vn,, Pl - € : 

with identical expressions for D"« v and D x r f 4 \. These rela- 
tions determine the earliest or latest time a conducting path 
from node i to Vdd can be established if at the beginning of 
the cycle no such path exists. Similarly, to determine 

D 7 a r ii V > D ) a *iu S/Vui we ^ P aths from the channel node 
(node t) to ground (node 0) with the target transistor off. In 
cases in which there are no paths, we set early arrival times 
to t cyci* and late arrival times to 0. 

5 Results and discussion 

We present static timing results from SOISTA for three ex- 
amples (of increasing complexity) and compare these with 
the results of circuit simulations in which vectors are chosen 
both to correctly sensitize the delay path in question and 
to match the assumed signal probabilities. At the 2.5-V 
supply, the accessibility results exactly match the detailed 
steady-state results because of the 6-state pinning; there- 
fore, separate accessibility results are not presented for the 
larger supply. 





Figure 3: Chain of identical inverters simulated with the 
waveform shown in the inset. 



5.1 Inverter chain 

The first is a chain of identical inverters as shown in Figure 
3 stimulated with the periodic waveform shown in the inset. 
The input waveform repeats every 10 nsec and is equivalent 
to a signal probability of 0,1 on the input of even-stage in- 
verters and 0.9 on the input of odd-stage inverters. Before 
t = 0, even (odd) stages have a zero (one) on their input. 
We consider results at two supplies, 2.5 V in Figure 4 and 
1.0 V in Figure 5. In Figure 4{a), we present the body volt- 
ages for the FETs as determined by circuit simulation for 
both even (in solid lines) and odd (in dashed lines) inverter 
stages. We notice that initially the even-stage FETs have 
higher body voltages than the odd-stage FETs. This is be- 
cause the even-stage FETs have their gates held low before 
t = 0 with the pFET (nFET) in state 2 (5), while the odd- 
stage FETs have their gates held high before t = 0 with 
the pFET (nFET) in state 5 (1). r b > r a for the nFET 
and ra > r$ for the pFET. As shown in the stage delay 
results in Figure 4(b) (even stage delay in solid, odd-stage 
delay in dashed), this gives the even stages initially a longer 
rise time but smaller fall time than the odd stages. Ini- 
tially, then the input pulses are stretched [16] as they move 
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down the inverter chain. As switching begins, state 6 be- 
comes accessible and the body voltages become pinned in 
about a microsecond (the time scale of re) to a&. We also 
note that, in this case, the steady-state nFET (pFET) body 
voltage is slightly less (more) than the initial-condition min- 
imum (maximum) value. 6-state pinning means that there 
is no difference between the steady-state body voltages of 
the FETs in even and odd stages. In steady-state, then, the 
rise and fall delays become the same in the even and odd 
stages and there is no pulse stretching. 



pfCT 




Figure 4: Inverter chain results at 2.5-V supply, (a) Body 
voltages from circuit simulation for even (solid lines) and 
odd (dashed lines) inverter stages, (b) Inverter delay. 



delays bound those determined by dynamic simulation. All 
of the uncertainty in this case is associated with the loading 
uncertainty, since the body voltage is nearly precisely pre- 
dicted by SOISTA. The delays determined by SOISTA with 
accessibility analysis are also noted in Figure 5(b). The 
initial-condition stage delays also bound the actual initial- 
condition stage delays of the dynamic simulation. In this 
case, the uncertainty comes from both the loading uncer- 
tainty and the body- voltage uncertainty. 




Tim* (mc) 



Figure 5: Inverter chain results at 1.0-V supply, (a) Body 
voltages from circuit simulation for even (solid lines) and 
odd (dashed lines) inverter stages, (b) Inverter delay. 



We now consider how these results compare with SOISTA. 
Early and late arrival times at the input of the chain, both 
rise and fall, are set to 100 psec. The rise and fall times 
at the input of the chain are also set to 100 psec. The di- 
amonds on the right vertical axes in Figures 4(a) and 5(a) 
correspond to the steady-state body- volt ages estimated by 
SOISTA for an input signal probability of 0.1, propagated 
as 0.9 to the inputs of the odd stages. State 1 and state 
2 for the pFET and nFET, respectively, are shown as solid 
diamonds. State 5 for the pFET and nFET are shown as 
hollow diamonds. The minimum and maximum values of 
the steady-state body voltage for these cases are indistin- 
guishably close because there is no path delay variation to 
produce any significant differences between the early and 
late arrival times. At 2.5 V supply, 6-state pinning dur- 
ing steady-state causes the detailed body voltages to be the 
same for even and odd stages, while at the 1.0-V supply, the 
absence of 6-state pinning causes disparate steady-state volt- 
ages between even and odd stages. The diamonds slightly to 
the left of the right vertical axis in Figure 5(a) are the min- 
imum and maximum body voltages determined by SOISTA 
from accessibility analysis. These numbers bound the full 
steady-state body-voltage variation of both even and odd 
stages. On the left vertical axis, we show the minimum and 
maximum initial-condition body voltage values estimated by 
SOISTA. These bound almost precisely the initial condition 
body voltages observed in the dynamic simulation, since in 
this case, the dynamic simulation covers all the accessible 
state states. In Figures 4(b), 5(b), the SOISTA steady-state 
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Figure 6: 16- bit ripple-carry adder circuit: (a) the associ- 
ated full-adder circuit and (b) the connection of the full- 
adder cells. The inset of (a) gives the waveform applied for 
dynamic simulation. 



5.2 Ripple-carry adder 

The next example we consider is a static 16- bit ripple carry 
adder shown in Figure 6(b). Figure 6(a) shows the compo- 
nent full-adder circuit. In the circuit simulation, we use the 
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input waveforms shown in the inset of Figure 6(a), which 
sensitizes the critical path of this circuit, the carry chain. 
The "A" waveform is applied to each A input and the a B" 
waveform is applied to each B input. The "C" waveform 
is applied to the Cxn input of the ripple-carry (see Figure 
6(b)). These waveforms correspond to signal probabilities 
of 0.5 on the A and B inputs of each full- adder cell and 0.5 
on the C input of each full- adder cell. For these input signal 
probabilities, the signal probability of Cout is 0.5, so that 
each cell sees identical switching activity. Before t — 0, A 
and C are high and B is low for each cell. Figure 7 shows 
the results for a supply voltage of 2.5 V. 
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Figure 7: Carry chain results at 2.5-V supply, (a) Body 
voltages for Ml and M2. (b) Full-adder cell delay. 

In Figure 7(a), we compare the body voltages of tran- 
sistors Ml and M2 of Figure 6(a) with the SOISTA initial 
condition and detailed steady-state results. (The detailed 
steady-state results match the accessibility steady- state re- 
sults because of the 6-state pinning.) The Ml curve labelled 
5 (2) and M2 curve labelled 1 (5) correspond to the case 
in which C is high (low) for a given cell. The steady-state 
values match almost exactly the values determined from cir- 
cuit simulation. They have negligibly small uncertainty be- 
cause the difference between the early and late arrival times 
is a small fraction of * C ycie. The initial-condition body- 
voltage values bound the circuit simulation results. The 
simulation results match exactly the lower bound for tran- 
sistor M2 because the simulation begins with M2 in state 
1, which determines the minimum possible initial-condition 
value. Ml, by contrast, begins in state 5 which is neither the 
minimum (state 3) nor maximum (state 4) accessible state 
(A»tatic = {1,3,5,4}). Therefore, the Ml initial condition 
values from simulation lie within the SOISTA-determined 
ranges. Figure 7(b) shows the complete stage delay of the 
full adder cell from C to Cout for both rising and falling 
Cout. The steady-state and initial-condition delays from 
SOISTA bound the simulation values. The uncertainty of 
the steady-state values derives entirely from the loading un- 
certainty. 

5.3 4-2 compressor 




Figure 8: 4-2 compressor circuit 
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Figure 9: Body voltages for Ml and M3 of 4-2 compressor 
circuit at 1.0- V supply. 



The last, and most complex, example presented here to 
compare the SOISTA and circuit simulation results is a 4-2 
compressor circuit from a tree multiplier design [17] (see Fig- 
ure 8). We specifically consider the long delay path, 14- A p- 
C-D-E-F-S, and the short delay path, I2-C-E-S. In Figure 9, 
we present body- voltage results for two FETs {Ml and M 3) 
along the long path at the 1.0- V supply. The dashed (solid) 
curves for MZ correspond to the case in which the gate is 
high (low) and the FETs are in state 5 (2). The dashed 
(solid) curves for Ml correspond to the case in which the 
gate is high(low) and the FET is in state 1 (5). The dia- 
monds on the left vertical axis give the SOISTA-determined 
initial condition body voltages for these same conditions. 
We note that at the dc state established in the circuit simula- 
tion before t = 0, II = 13 = 0, 12 = 14 = Cin = 1.0 V. At these 
values, we find C = 0, B =: 0. Therefore, M3 is in state 4 and 
Ml is in state 1. For the static accessible sets of Ml and 
M3 {Astatic = {5,1,3,4} for M3 and Ataxic = {1,2,4,5} 
for Ml), these represent the states of highest and lowest ref- 
erence body voltage, respectively. Therefore, these curves 
match the maximum (for M3) and minimum (for Ml) body 
voltage values calculated by SOISTA. On the right vertical 
axis, we show the steady- state detailed body voltage associ- 
ated with the longest path sensitization. As expected for the 
late-mode case, these are slightly less than the simulation 
values for the nFET and slightly greater than the simula- 
tion values for the pFET. We also show the body- volt age 
estimates from accessibility analysis, which quite expect- 
edly bound the detailed results. In Table 1, we compare 
the initial-condition delays determined by circuit simulation 
with the SOISTA-determined initial-condition delays. Table 
2 does a similar comparison between the steady-state delays 



537 







SOiSTA 


initial 




SPICE initial 




Rise | 


Fall 


short-path f 


long-path 




max 


nun 


max 


nun 


rise 


fall 


rise 


fall 




100 


100 


100 


100 


100 


100 


100 


100 


14 


too 


100 


100 


100 


100 


100 


100 


100 


Ap 


285 


229 


398 


284 






235 


315 


C 


645 


318 


576 


320 


326 


379 


476 


432 


D 


878 


450 


946 


501 






723 


553 


E 


1341 


531 


1200 


537 


629 


762 


785 


943 


F 


1602 


661 


1366 


618 






1112 


936 


S 


1842 


681 


1905 


662 


1064 


837 


1104 


1174 



Table 1: SOISTA-determined initial-condition delays versus 
circuit simultaion delays. 
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Table 2: SOISTA-determined detailed steady-state delays 
versus circuit simulation delays. 



determined by circuit simulation (after more than 50,000 
cycles of simulation) and the SOISTA-determined detailed 
steady-state delays. In all cases, the SOISTA-determined 
delays bound the SPICE delays. One should also notice the 
considerable reduction in uncertainty between the initial- 
condition and detailed steady-state delays, as the compo- 
nent of this uncertainty due to body voltage variation is 
noticeably reduced. (The remaining uncertainty is due pri- 
marily to delay- path variation and loading uncertainty.) 

6 Conclusions and future work 

In this paper, we have presented a circuit-focussed model of 
the floating-body potential of PD-SOI FETs. This model 
allows one to determine the body voltage and its associated 
uncertainty, depending on knowledge of the switching ac- 
tivity of the FETs in question. Four types of estimation 
are possible depending on switching assumptions and the 
amount of information known about the logical and tempo- 
ral environment of the circuit under analysis. We have in- 
corporated this model into a prototype transistor-level static 
timing analysis engine to demonstrate the impact reduced 
body- volt age uncertainty can have on performance evalu- 
ation. We find that the body- voltage uncertainty can be 
significantly reduced with fairly conservative assumptions 
about switching behavior. 

Future work will include incorporating these body volt- 
age estimates into transistor-level static noise analysis. In 
addition, we intend to consider design techniques whereby a 
normally inactive block could be periodically stimulated to 
keep it "primed" so that when it is eventually exercised, it 
has more tightly predictable body voltage variation. This is 
similar to some of the circuit techniques which attempt to 
force discharge of the body during "non-critical" periods of 
circuit operation (e.g. precharge in dynamic logic) to reduce 



parasitic bipolar leakage. In many ways, this could also be 
viewed as analogous to DRAM refresh. More work will be 
required to determine the necessary frequency and nature of 
this pattern. 
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